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Abstract 

We propose to study the superconducting pairing of quarks with the formation of 
the diquarks as well as the quark-antiquark pairing in QCD by means of the functional 
integral technique. The dynamical equations for the superconducting order parameters 
are the nonlinear integral equations for the composite quantum fields describing the 
quark-quark or quark-antiquark systems. These composite fields are the bi-local fields if 
the pairing is generated by the gluon exchange while for the instanton induced pairing 
interactions they are the local ones. The expressions of the free energy densities are 
derived. The binding of three quarks is also discussed. 

The superconducting pairing of quarks due to the gluon exchange in QCD with the formation of 
the diquark condensate was proposed by Barrois' 1 ' and FrautstriJ 2 ! since more than two decades and 
then studied by Bailin and Love' 3 !, Donoglue and Sateesh", Iwasaki and Iwado' 5 '. Recently, in a 
series of papers by Alford, Rajagopal and Wilczek' 6 ] , Schafer and Wilczek^l , Rapp, Schafer, Shuryak 
and VelkovskyI 8 ! , Evans, Hsu and Schwetz' 9 !, Son' 10 !, Carter and Diakonovl 11 ! and others' 12-21 ' there 
arose a new interest to the existence of the diquark Bose condensate in the QCD dense matter-the 
color superconductivity. The connection between the color superconductivity and the chiral phase 
transition in QCD was studied by Berges and Rajagopal' 22 ' , Harada and Shibata' 23 ! . There exists 
also the spontaneous parity violation, as it was shown by Pisarski and Rischke' 24 ' . 

For the study of many-body systems of relativistic particles with the internal degrees of freedom 
as well as with the virtual creation and annihilation of the particle-antiparticle pairs the functional 
integral technique is a powerful mathematical tool. This method was applied to the study of the 
color superconductivity as well as the quark-antiquark pairing in a series of papers of the group of 
authors in Hanoi' 25-28 ! . In my talk I present a review of these works. 

Denote ip A the quark field, where A — {a, a, i) is the set consisting of the Dirac spinor index a — 
1, 2, 3, 4 and the color and flavor indices a = 1, 2, 3...N C and i = 1,2, 3, ....Nf. The internal symmetry 
groups are assumed to be SU (N c ) c and SU (Nf)j.We work in the imaginary time formalism and 
denote x — (t, x) a vector in the Euclidean four-dimensional space-time. The partition function of 
the system of free quarks and antiquarks with the chemical potential [i and at the temperature T 
can be expressed in the form of the functional integral 

Z = J [Di>] [D$] exp | - j dxi, A (a?) D^ B (x) J (1) 
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with brief notations 



P 

Jdx^JdrJdx, /?=£f. D A= 7 4 (^:-m) +7V + M 



(2) 



where k is the Boltzmann constant and M is the bare quark mass. We start from the study of the 
system with the direct instanton induced four-fermion couplings and then consider the case of the 
pairing due to the gluon exchange. 

We write the expression of the interaction Lagrangian of the direct four-fermion couplings in two 
different forms 

Li n t = 2^ {x)tp {x)Vca V , d(x)iP ' B (x) ( 3 ) 

and 

L int = \^ A {x)^ B {x)U^{x)^ D (x) (4) 

convenient for the study of the quark-quark or quark-antiquark pairing, resp. The notations 
V^a an d ^ac are related each to other 



v bd _ jtBD 
V CA — u AC ■ 

These coupling constants are antisymmetric under the interchanges of the indices. 

In order to study the quark-quark pairing we use the interaction Lagrangian in the form (3) and 
have the partition function 

Z = J [Dip] [D$\ exp | - J dx$ A (x)D B ^ B (x) J 

.expji J dx^ix^ix^f^Dix^Bix)^ . (5) 

AC 

Introduce the composite local bi-spinor fields ^db(x), $ (x) describing the diquarks and consider 
the functional integral 

Z*'* = J [D<$>][D$] exp j-1 J dx^'ix^f^DBix)^ . (6) 

Applying the corresponding Hubbard-Stratonovich transformation, we rewrite Z in the form of the 
functional integral over these composite local fields 



Z=^= /[£>*][!>*] exp {S cS [*,$]}. (7) 

The field equations 



2 



8S og [$, g 
5^ AC {x) 



= 



(8) 



are the dynamical equations for the quark-quark pairing. They can be written in the following 
manner 

A DlBl On) = V£l°l | J dx 2 S% ( Xl - x 2 )A C2 a 2 (x 2 )Sil ( Xl - x 2 ) 

- j ' dx 2 J dx 3 J d Xi Scl{xi - x 2 )A C2 A 2 (x 2 )S^l(x 3 - x 2 ) 



-A 3 C 3 



A CA (x) = Vgf* DB (x), A" u {x) = Vg2* AC (x). 



(9) 



(10) 



•A" 3 " 3 (x 3 ) S% (x 3 - x 4 ) A CiAi (xi) SH ( Xl - x 4 ) + } 
where S%(x — y) are the propagators of the free quarks and 

— BD , 

Consider the special class of the constant solutions 

a ir BD 
Ac a — const, A — const. 

Denote 5^(p, e m ) the Fourier transforms of the free quark propagators, e m — |(2m+l), and 

introduce A, A and S(p,e m ) the matrices with the elements Ac a, A B ° and S'f (p,e m ). Then we 
have following matrix equation for the quark-quark pairing 



DB 



V AC 
V DB 



5? (£)*/* 

1 



l + S(p,e m )AS T (-p,-e m )A 
The corresponding free energy density equals 

F = \ E 7^3 / d ?\ Tr p(p,e m )A5 T (-p, -e m )± 



-S(p, e m )AS T (-p, -e m ) 



(11) 



CA 



da- 



l + S(p,e m )AS T {-p,-e m )A J 1 + aS(p, e m )AS T (-p, -e m )A 



(12) 



In order to study the quark-antiquark pairing we use the interaction Lagrangian in the form (4) 
and start from the partition function 



. exp 



exp 



—A 



dx^ A (x)D^ B (x)^ 



dxiP {x)^ B {x)U^4, {x)i, D {x)\. 



(13) 
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Introducing the hermitian composite fields &b(x) as well as the functional integral 



(14) 



and applying the corresponding Hubbard-Stratonovich transformation, we rewrite Z in the form of 
the functional integral over these new composite fields 

Z= J| |[£$]exp{S cff [$]}. (15) 

The field equations 

are the dynamical equations for the quark-antiquark pairing. They have the explicit form 

&S{x) = Ug%Gi(x,x) (17) 

where 

AS(x) = Ug%*i(x) (18) 
and G A (x, y) are determined by the Schwinger-Dyson equation 



G B A (x, y) = Si(x-y)- f dzS c A (x - z) A B (z) G B D (z, y) . 



(19) 



It is easy to verify that G A (x, y) are the two-point Green functions of the quarks in the presence of 
the quark-quark pairing. In the case of the constant order parameters 

= const 

the Green functions G A (x, y) depend only on the coordinate difference x — y. Their Fourier trans- 
forms G A (p,e m ) are determined by the algebraic equation 

G B A (p, e m ) = S% (p, e m ) - S C A (p, e m ) AgGg (p, e m ) . (20) 

Introducing the matrices G(p,e m )and A with the elements G A (p,e m ) and A^, we obtain the 
matrix equation 

1 +A. (21) 



G(p,e m ) S(p,e m ) 
The free energy density of the system equals 



(22) 



i 

1 



da- 



2l + A(p,e m )S(p,e m ) J 1 + aA(p,e m )S(p,e m ) 



o 



4 



For the study of the quark-quark pairing due to the gluon exchange we start from the partition 
function in the form 



Z = 



J [£ty][Ity]exp dxil) A {x)D*i> B {x) 
.exp^jdxj dy^ A {x)^ C (y)V^ (x - y) iI>d(vWb(*)} > ( 2 3) 



where 



(z-y) 



c6 1 rfe rrf 



(24) 



A/ are the Gell-Mann matrices of the color symmetry group. In order to describe the diquark 

AB 

systems we introduce some bi- local composite fields $cr> {y, x) , $ (x, y) . Consider the functional 
integral 

Zp = y^^expj-1 jdxj dy^ AC {x,y)Vi^{x-y)^ DB {y lX )Y (25) 

Applying the corresponding Hubbard-Stratonovich transformation, we rewrite the partition function 
in the form 



Z = / [ D *\ l D ^} ex P i S M®, *] } • (26) 

z o J 

The field equations 

(5$ (a;, j/) 

are the dynamical equations for the quark-quark pairing. The bi-local composite fields §db(v — 
\t • 

x), <& {x — y) depending only on the coordinate difference, or the linear combinations 

Aca (y-x) = VSa (x - V) *D B (y - x), A B ° (x-y) = v£f (x - y) - y), (28) 

„ ^_BD 

are the order parameters of the superconducting phase transition. Denote Ac a (p, s m ) and A (p, e m ) 

their Fourier transforms and A (p, e m ),A (p, e m ) the matrices with these elements. Then the field 
equations (27) become the extended BCS equations for the color superconductivity induced by the 
gluon exchange 
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A C A (P, e m ) = -g 2 (7„ ® ( 7/1 ® A,)£ .- £ ^3 j dq 



(e n - £ m ) 2 + (p - q) 2 



1 + S (q, e n ) A (q, e„) 5" T (-q, -e n ) A (q, e n ) 



■S (q, £„) A (q, e n ) S T (-q, -£„) 
The free energy density of the system equals 



(29) 



F[A,A] 



-T — 



dpTr S (p, e m ) A (p, e m ) S T (-p, -e m ) A (p, e„ 



1 



1 + S (p, e m ) A (p, e m ) S T (-p, -e m ) A (p, e m ) 
1 

1 



(30) 



— / da- 



1 + aS (p, e m ) A (p, e m ) S T (-p, -e m ) A (p, e m ) 
In order to study the quark-antiquark pairing we use the partition function (23) in another form 

Z = y[£»^[£»^]exp|- j ' dx$ A {x)D^ B {x)} 

.expjl Jdxj dy^ A (x)^ B (y)U^ (x - y)^' (y)^ D (x)Y (31) 

U%g(x-y) = -Vg'f(x-y). 

The quark-antiquark systems are described by the bi- local hermitian composite fields $b(x, y). 
Introducing the functional integral 



Zt = f [£$]cxp j-i Jdxj dy*i(x,y)U%g(x-v)*%{v,x)} (32) 
and applying the corresponding Hubbard-Stratonovich transformation, we rewrite Z in the form 

Z = §k |[^]exp{5 cff [$]}. 



The field equations 

(5$g (y,x) 

are the dynamical equations for the quark-antiquark pairing 



(33) 



(34) 
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*%{v,x) = UB%{x-v)Gi{y,x), (35) 
where Gg (y, x) are determined by the Schwinger-Dyson equations 

Gi (y, x) = Si(y-x)- J dx' J dy'sg (y - y') A% (y' , x 1 ) G% (x' - x) (36) 

and 

A B A {x,y) = U%g{x-y)$ c D {y,x). (37) 

(y, x) are the two point Green functions of quarks in the presence of the quark-antiquark pair- 
ing. In the special class of the functions A^ (x — y) depending only on the coordinate difference these 
functions are the order parameters of the system. Denote A^ (p,£ m ) their Fourier transforms and 
A(p, e m ) the matrix with these elements. The dynamical equations (35) are 



>c 



Ag(p,e m ) = 5 2 ^(7 M ®A / )JJ(7 Al ®A / ); 

I If 1 ~ 

-rY^TTTTs / d % 7 -2G^(q,e„), 

P ^ (2tt) J (e n - e m ) + (p - q) 



(38) 



and the Schwinger-Dyson equations become 

1 1 



G(q,e m ) S(q,s m ) 
The free energy density of the system equals 



+ A (q, e m ) . (39) 



F t A l = \ E ^3 / d P Tr [ A (P- Cm) S (p, e m ) 



— f da— 

2l + A(p,e„ l )S'(p,£ m ) J 1 



+ aA (p,e m )S (p,e m ) 



(40) 



The existence of the non- vanishing order parameters Ac A: A or Ac a (y — x), A (x — y) 
which cannot be the singlets of the color symmetry group would mean the spontaneous breaking of 
the color symmetry! 29 ] . 

The physical vacuum in QCD is the limiting case of the system with vanishing density (p = 0) 

at vanishing temperature (T = 0). In this limit the order parameters Ape, A and Ag are the 
vacuum expectation values of the corresponding quantum composite fields - the new Higgs fields. 

The functional integral method was applied also to the study of the binding of three quarks due 
to some effective six-fermion (non-local, in general) interactions of quarks. The systems of dynamical 

Aj^C 

equations for the composite fermionic quantum triquark fields ^ abc an d VP were derived. The 
constant solutions of these equations in the case of the direct six-fermion couplings of quarks or 
their solutions depending only on the coordinate difference in the case of the non-local effective 
six-fermion interactions are the anticommuting order parameters. 
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